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A PARTICULAR CLASS OF SOLUTIONS OF A SYSTEM OF 
EIKONAL EQUATIONS 

THIERRY CHAMPION AND GISELLA CROCE 

Abstract. In this paper wc study the foUowing Dirichlot problem: 

du 

1 a.e. inQ, i = 1, . . . , N 

on an, 



dxi 



rf\ , where O, is an open bounded subset of R^ and u S VF^'°°(n). This problem 

has infinitely many solutions; our aim is to isolate those which minimize some 



functional involving the discontinuity sets of ^^, for i = 1, . . . ,N . 



1. Introduction 

In this paper we are interested in the foUowing system of eikonal equations 
du 



dxi 



— 1 a.e. inil, i = 1, . . . , A^ 



(1.1) 
^ \ [ M = on do., 

C^ . where fi is a connected open bounded subset of R^ with Lipschitz boundary and 

Cn ' u S W^'°°(ri) n Co{Vl). This problem admits infinitely many solutions (see [5] and 

[S] for example for a proof). An interesting question is then to isolate a particular 
lO ' class of solutions which have some additional properties, since such an analysis 

^? , could be useful to select and characterize a prefered solution. 

There exist some evident difficulties to this kind of question. For example the 

lack of convexity of the set &{yi) of the solutions of p.ip implies that the natural 

functionals 

>< ■ " -> / 1^1" P ^ 1' 

H ; . . . .^ . . 

5^ ! have in general neither a minimizer nor a maximizer over &{Vt). Indeed any min- 

imizing sequence converges to 0, which does not belong to ©(^2), and the cluster 
points of maximizing sequences are the distance function to dVl (in the l^ norm of 
M^) or its opposite, which usually do not verify (|l.ip (see remark [5^ . 

Another natural way to select a particular solution in &{Vt) would be to use the 
notion of viscosity solution (for the definition and further details see [5] ) . Indeed 
one could study the existence of iy^'°°(ri) viscosity solutions of a problem like 

f FiDu) = a.e. inf7 
\ u — VI onoil 



1991 Mathematics Subject Classification. 28A75, 34A60, 35F30, 35A15, 49J45. 
Key words and phrases, almost everywhere solutions, eikonal equation, functions of bounded 
variation, direct methods of the calculus of variations. 

The research of Gisella Croce was partially supported by the Fonds National Suisse. 

1 



2 THIERRY CHAMPION AND GISELLA CROCE 

where F : M.^ — > M is such that F{x) = if and only ii x — {xi, . . . , xn) verifies 
\xi\ = I y i = 1,..., Xn- This is indeed possible in dimension A^ = 1, where the 
only viscosity solution of the problem 

, ,-. , - . a.e. inl — 1, If 

(1-3) <! :./.. ^... _^ ^ ' ^ 

is u{x) = 1 — |a;|. Nevertheless when N >2 there is in general no viscosity solution 
of ()1.2p . unless fl is convex and the normals to dil are in the set E — {x : \xi\ = 
1/\/N}] we refer to [S] and [13] for further details. 

In the literature we find few articles with the purpose of isolating particular 
solutions of problem (|l.ip and selecting one. In 7 and [5] the authors study 
numerically a variational problem over the set of non- negative solutions of (jl.ip in 
the 2-dimensional case: they obtain a maximizing sequence for the problem 




(1.4) 



sup \ u, M > 0, u e 6{rt) > 



through the numerical minimization of the functional 






Mi~>- / u+- I \Vur + - I |A 



2 / ' ' 2 / ' ' 2e 



N 



= 1 



d 



u 



dxi 



- 1 



Unfortunately, as we said above, there is in general no optimal solution for the 
variational problem (|1.4|) . since the limit for any maximizing sequence is the distance 
to dil (in the l^ norm). 

A different approach is followed in |10j where the authors construct an explicit 
solution. Indeed they define a special partition of $7 made up of convex sets whose 
normals to dO, are in E\ over each of these sets they consider the relative viscosity 
solution. In this way the constructed solution has local properties (which depend 
on the chosen partition), but no global ones. In any case we want to point out that 
this study is made in a more general framework than problem (jl.ip . 

Our purpose in this work is to characterize a particular class of solutions of 
(|l.ip through a variational method, i.e., as the optimal solutions of a variational 
problem. We shall define over &{Vt) a functional involving the discontinuity sets 
of the gradient of the solutions of (|1.1|) and we minimize it. Let us briefly explain 
the kind of results that we shall present in this article. We distinguish two cases, 
according to the properties of 57. In the case where 917 is composed by a finite set 
of faces with normals in E (see section [3]), we consider the functional 

N 

i=l 

(where Ja^ is the jump set of 4^) and we show that it admits a minimizer. 

In the general case, that is, when J7 is any connected open subset of R^ with 
Lipschitz boundary, we cannot minimize the discontinuity sets of Du for u G 6(57) 
in the same way. Indeed the previous functional is in general infinite over ©(17), 
as shown in section H) Consequently we construct a positive "weight" function 
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h e Co{Q) such that the functional 



N 

v-^J2 h{x)d'H^-\x) 



'Jo. 

admits a minimizer over 6(51). 

We also study the relationship between the viscosity solution (when it exists) 
and the optimal solutions of the minimization problem: we show that they coincide 
if iV = 1, 2 (see Proposition EH). 

Our study allows us to isolate particular solutions of problem (jl.ip , but unfor- 
tunately we are not able to say if it selects a particular one, in the sense that the 
variational problems we propose could admit more than one solution. In section [5] 
we then present a method to select a function among the optimal solutions of the 
variational problems that we studied (see 4.). 

Finally we refer to a recent work ^ which does not address the problem of the 
selection of a particular solution of (|l.ip but is also concerned with some local 
properties of the discontinuity sets of Du, u E &{il). 

2. Notations 

Throughout this paper J7 denotes a bounded connected open subset of M^ with 
Lipschitz boundary. We will denote by {xi)i<i<N the coordinates of a point x e R^, 
and (x, y) the usual scalar product of x, y G M.^ . The norm ||a:|l;i of x G R^ stands 
for the usual l^ norm, that is ||a;||;i = J2i=i l^d ■ 

In the following, we shah use the short notation Wo'°°(rj) for W^'°°{n)nCo{n). 
When considering a function v € VFQ'°°(r2), we shall always assume that we are 
handling its continuous representative; notice that, due to the hypothesis on J7, 
such a function w is a Lipschitz continuous function. 

Let us briefly recall some properties and notations of the spaces BV and SBV for 
future use (in any case we refer to and to |Tl] for further details). If w € BV{il.), 
we write HwH^y/j^N for the usual BV norm of w defined by 

ll"'llsv(o) = IMl^u) + \Dw\{n) . 
Moreover we will use the usual weak* convergence in BV defined as follows (see 

my- 

Definition 2.1 (weak* convergence in BV). Let (m„)„,m g BV{fl). We say that 
{un)n weakly* converges to u in BV{il) if u„ — > u in L^{Q,) and the measures Dun 
weakly* converge to the measure Du in A^(il,M^), that is 

lim / ipdDun — / ipdDu V (p E Co{^) ■ 

A BV function belongs to SBV when the Cantor part of its derivative is zero. 
This implies that the measure Dw of an SBV{^) function can be decomposed in 
the following way: 

Dw = D^w + D^w 

= \/wC^ + {w+ ~w-)i^^'H^-\j^ = \/wC^ + [w]i^^'H^~\j^ 

where D'^w is the absolute continuous part of Dw with respect to the Lebesgue 
measure £^ in M^, D^w = [wjVwT-i^ ~^ [Jw is the jump part of Dw, 'H^~^ the 
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Hausdorff measure of dimension iV — 1, w+ and w^ denote the upper and lower 
approximate limits of w, Jw the jump set of w and Vw its generalized normal. 

Referring to our problem, we will denote by S(J7) the set of W^'°°{ft) solutions 
of dUI]), that is 



&{n):=\vew^-^{n) -. 



dv 



dxt 



1 a.e. in Q,i ^ 1, . . . , N 



In this paper we will consider the piecewise C^ elements v of ©(fi) and more 
generally those elements v of &{fl) such that Dv S SBViod^)'^ ■ To define a 
piecewise C^ function, we first recall the definition of Caccioppoli partitions and 
piecewise constant function (we refer to [5] for further details): 

Definition 2.2 (Caccioppoli partition). Let I C N. A partition {Ei}i^i of il is a 
Caccioppoli partition if \^ \DxEi\{^) < co- 

In the above definition, XA denotes the characteristic function of a set A, that 
is xa{x) = 1 when x £ A, xa{x) = when x ^ A. 

Definition 2.3 (Piecewise constant function). We say that u : fi -^ R is piecewise 
constant in fi if there exist a Caccioppoli partition {_Ei}ig/ of D, and a function 
i : / ^- K such that u ~ >, iiXEi ■ 

Notice that a piecewise constant function is an SBV function (see [2] for further 
details) . 

We define the set of piecewise C^ functions by 

Cp^{^) := Iv e W^'°^{n) : — ^ is piecewise constant in n,i ^ I, . . . ,n[ . 

As has already been said, we shall handle the solutions of (ll.l[) such that Dv G 
SBVioci^)^ ■ We now want to write their distributional gradient: for any index 
i € {1, . . . , N} we use the short notation J^ instead of J s^ . Therefore for a function 

V e6{n) such that Dv e SBViod^)^ one has 



(2.1) ^[b~] = ^''i^'^ L-^*" ^"^^ Vi=l,...,iV 

for any open subset w C ZU C il. 

In the rest of the paper, the distance in l^ norm to a subset of R^ will play an 
important role: if A is an open bounded subset of M^, we define dji {■ ,A) : R^ — > R 

by 

VxeM^, diiix,A) := min{||a; - y||,i : y e A} . 

3. A PARTICULAR CASE 

In this section we will assume that fl satisfies the following property: 

the boundary dfl is the union of a finite family of polyhedrons 
(H) ^ of dimension A'^ — 1, each of which is included in a hyperplane 

with normals in E, 
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where E := {x € R^ : Vi, 1x^1 = l/\/iV}, and we will minimize the H^ ^ measure 
of the discontinuity sets of ^^, for v G &{n). We define 

We now state the main result of this section. 



Theorem 3.1. Let il be an open bounded connected subset o/ R^ which satisfies 
hypothesis (H). Then the variational problem 



N 

V^TjAT-l/ 



(Pi) mf <^ T{v) := }^H^-\J^) : v e Sp^i^) V . 

has finite value and admits an optimal solution. 

Remark 3.2. Under assumption (H), problem (Pi) always has at least two optimal 
solutions: indeed, if u is a solution then —u is also a solution. 

Remark 3.3. When ft fails to satisfy hypothesis (H), then problem (Pi) is not well 
posed as we will show in section |4l 



In the proof of Theorem l3.1l we will use the following lemmas. 

Lemma 3.4. Let fl be an open bounded connected subset of M.^ with Lipschitz 
boundary. Then 

-dii{-,dn) < V < dii{-,dn) onTi 
for every function v £ &{il). 

Proof. Let v G 6(il). We denote by A a set of zero Lebesgue measure such that v 
is differentiable on il \ A with 



dv 



= 1 for all i. 

Take a: G fi, and let y g di^ be such that d;i {x, dil) ~ \\x — y||;i . For any t E ]0, 1] 
we define the map jt '■ [0, 1] — > fi by 7t(s) :— sx + {l — s)x^, where x* := tx + {l — t)y. 
For a fixed t s]0, 1[ and e > 0, we apply Lemma 3.2 in |4l to get the existence of 
a Lipschitz continuous path 74^^ : [0, 1] — >■ fi \ A such that 

7t,e(0) = x\ 7t,e(l) = X and |J7^^^ - 7t||Loo(o,i) < e. 

Then the function v o ^t.e is Lipschitz continuous on [0, 1] and we can compute 

v{x)-v{x*)^ {v o -ft^^y {s)ds ^ Vi;(7t,e(s)).74„(s)ds. 
Jo Jo 

Now ^^(7t,s(s)) = 1 for all i and so 

v{x)-vix*)< [ \\-f'tJs)\\ ds<Ne+ [ His)\\i,ds = Ns+\\x-x'\\^,. 
Jo Jo 

Letting t and e go to this yields v{x) < dii{x,dQ). In an analogous way one can 

prove that —v{x) < rf/i (x, 957) D 

Lemma 3.5. Let il be an open bounded connected subset of M.^ which satisfies 
hypothesis [H). Then (i;i(- ,(9^2)) belongs to £pw[Q) and J^(d/i(- ,9^2)) < +00. 

Remark 3.6. Notice that for a general bounded open subset i7, d/i ( • , 957) does not 



necessarily satisfy 



dv 



1 a.e. in J7 for all i. Indeed, it is sufficient to consider 



dxi 

ri — (0,1)^ C M^ : we have that dii{{xi,X2),d^) = X2 in the set T = {(xi,a;2) 
<xi < 1/2, < a;2 <xi}. 
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Proof. Since J7 satisfies hypothesis (i/), we can write 951 = M i^j where each face 

Fj is a polyhedron of dimension iV — 1 included in a hyperplane Hj with normal 
in E. We shall assume that two faces Fj and Fk that intersect are not included in 
the same hyperplane of M^. 

Let i G {1, . . . , m}. We first notice that since the normal of the hyperplane Hj 



is included in E, one has 



-dA-.H^) 



— 1 a.e. in J7 for all i. Let us prove that 
dii{-,dVL) = min dii{-,Fj) also satisfies this system of eikonal equations in fl. 

l<j<rn 

We fix X in 51 and distinguish three cases: 

Case 1. Assume that d;i(x,9r2) = dii(x,Fj) for a unique index j. Then 
dii{x,Fj) is attained in the relative interior of Fj, so that in a neighborhood of 
X one has dii (• , dVt) = d;i ( • , Hj), and thus in this neighborhood d;i (• , dVl) is reg- 
ular and satisfies the system of eikonal equations (ll.ip . 

Case 2. Assume that dii {x, 951) = d/i {x, Fj) = d;i {x, F^) is attained in Fj n Fj. 
with j 7^ k. Then dii{x,Fj) = dii{x,Fi^) is attained in the intersection i/^ n TJ^,. 
We notice that the set of the elements y of M^ such that dji (x, Fj) = dji (a;, Ffe) is 
attained in Fj D Fk is included in the hyperplane Aj^k generated by Hj D Hk and 
£,j + ^k, where ^/ denotes a normal vector to Hi for ^ e {j, A:}. Since £''^(^j^/c) = 
this set is negligeble in view of the system of eikonal equations (11.11) . 

Case 3. Assume that d/i {x, 951) = d;i {x, Fj) = d/i (x, Fk) for j ^ k such that the 
hyperplanes ffj and i?fc are parallel and that case 2 is excluded. Then d;i (x, F) is 
attained in the relative interior of F for I S {j, fc}, and thus one has d;i (x, i/j) ~ 
dii{x,Hk)- As a consequence, x belongs to the set Bj^k '■= {y ■ dii{y,Hj) = 
dii{y,Hk)}, which is an hyperplane parallel to Hj and Hk, so that C^[Bj,k) = 0. 
This means that this set is negligeble in view of the system of eikonal equations 

dni). 

As a consequence of the above discussion, almost every x in 51 enters case 1, so 
that d;i ( • , 951) is a solution of ()1.1|) . Moreover, if z e {1, . . . , N^ then one has 

c |J[5in(Aj-fcUBj,fc)]. 

Since 51 is bounded and the sets Aj,k and Bj^k are hyperplanes of R^, we infer that 
the last union has finite H^~^ measure. This being true for all index i, we get that 
J"(d/i(-, 951)) is finite. D 

We can now turn to the proof of Theorem 13.11 

Proof of Theorem \3.1\ As a direct consequence of Lemma 13. 5[ we get that the in- 
fimum inf(Fi) is finite. We now apply the direct methods of the Calculus of Vari- 
ations to prove the existence of an optimal solution of (Pi). 

Let (u")„gN C £'ptu(51) be a minimizing sequence for (Fi). We deduce from 
Lemma [3.41 that (w")n is uniformly bounded in F°°(51). Moreover this sequence is 
uniformly Lipschitz continuous on f2 since 9f2 is Lipschitz and |Vu"| < vN a.e. 
in 51 for any n. Applying Ascoli-Arzela and Banach-Alaoglu Theorems, and up to 
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a subsequence, we can assume that u" — >■ v°° in Co{n,) and v^ — >■ v°° weakly* in 

W^'°°{n) for some w°° G W^'°°{n). 

Step 1: we show that v°° G £p^(f2). We first prove that ^^^ = 1 a.e. in J7. Let 

us fix z G {1, . . . , TV}. We observe that since (w")ra is a minimizing sequence, there 
exists a constant C > such that J^(w") < C for every n S N. This implies that 



(3.1) 



(9a;i 



= C^{n) + 2H^-^{jf) < £^{n) + 2C. 



BVin) 



Applying Theorem 16.11 to the sequence ^^ we have that ^^ — ?> gi weakly* in 



BV{fl) for some function gi e BV{Q,) (up to a subsequence). Since we already 



Moreover 



know that ^ ^ ^ weakly* in L°°{n), we infer that ^ = 9', 
f^ ^ ^ in LV17) so that 1^1 = 1 a.e. in n. 

OXi OXi ^ ' I OXi I 

To prove that v°° belongs to Cp^(51) it is sufficient to apply Theorem 16.31 to 



the sequence 



dv^ 



Indeed, 



dv" 



is a piecewise constant function for every n and 



dv^ 



dxi 



L°°(n) 



H^ ^(J" ) is uniformly bounded since F{v^) < C. As a conse- 



dv^ 



quence, up to a subsequence, ^j|- converges in measure to a piecewise constant 



function which is necessarily ^j— ■ This proves that v°° G £pw{^)- 

Step 2: we show that v°° is a minimizer of J^. Since {vn)n is a minimizing sequence 

for J^ it is sufficient to prove that 



(3.2) 



liminf H 



N-l 



{j: ) > n^'-'ij: ) 



for alii ^1,...,N. Thanks to (12.11) one has for any n G N U {+00} 



dx 

dv^ , dv 



in) = 2n^'\jf). 



Moreover we know that ^j- -^ ^^ weakly* in BV{il). Applying Theorem 
with / = 1 we get (|3?2|) . 



D 

We now study the relationship between the viscosity solutions of problem (jl.ip 
and the optimal solutions of problem (Pi). We restrict ourselves to the one and 
two-dimensional cases, where the geometric arguments are intuitive. 

In the one-dimensional case, it is evident that the viscosity solution u(x) = 1— |a;| 
and its opposite minimize the number of jumps of the solutions of (jl.3|) . as their 
derivatives have just one jump. For the 2-dimensional case, we recall the following 
result about the existence of a viscosity solution (see Example 1.3 and Theorem 3.3 
of [S]) . 

Proposition 3.7. Let F : M? ^ M. be a continuous function such that ^(^1,^2) = 
if and only if \^i\ = 1, i = 1, 2 and F{^i, ^2) < if \^i\ < 1 for i = 1,2. Let fl be an 
open bounded connected subset ofM.'^. Then there is no W^'°°{fl) viscosity solution 
of problem 

j F{Dv) =0 a.e. inn 

] V — on dfl, 

unless il is a rectangle whose normals are in the set E = {x — {xi, X2) G M^ : \xi\ = 
-7^ ,i = 1,2} . In this case the only viscosity solution is dii (• , dfl). 
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It also happens that d/i(- ,dfl) and its opposite are the only optimal solutions 
of (Pi) if 17 is a rectangle with normals in E, as the following result shows. Notice 
that if il C M^ is convex and satisfies property (H), then il is necessarily such a 
rectangle. 

Proposition 3.8. Let fl be a rectangle with normals in E. Then ±d/i(- ,951) are 
the only optimal solutions of (Pi). 

We will use the following lemma in the proof: 

Lemme 3.9. // P : M" — > K. is the projection (xi,X2) — >■ X2, and E is a measurable 
set o/K^ then 

[ H^{E n p-^{y})dy < H^{E). 

Proof. It is sufficient to apply Theorem 2.10.25 of [12^ with f = P, X = R'^,Y = R, 
k = and m = 1. D 



We are now able to prove Proposition lS.SI 

Proof. To simplify the notations, we set d := (i;i(- ,951). Let v £ £^^,(51) different 
from d and —d: we want to prove that 

2 2 

i=l 1=1 

We first claim that rO-{Jl) > H^iJ^)- To this end, we define 

^ := {f e M : 3 .s e M such that x := (s, t) e ft and -d{x) < v{x) < d{x)} . 
Since the functions v, d and ~d are continuous, the set A can be written as M Ji^, t,^[ 

where / C N is non empty, and tt < t[ for all i e /. We notice that 

H\Jl n (M X T)) = H^Jt n (M X T)) 

where l"" = R \ l^. We also notice that U^iJf n (M x dA)) = 0, so that in fact 

n^{Jl n (R X A")) > n\Jt n (r x A")). 

For L^-dlTnosX every t £ A, function w : s i— > v{s,t) is differentiable on ]a,6[ with 
derivative ±1, where a and b are such that ](a, i), (6, i)[= 51 n (R x {i}). Moreover, 
— d((s,t)) < w{s) < d{{s,t)) for some s s]a, &[, so that w' — -^ has at least two 
jumps in ]a, 6[. As a consequence, applying Lemma 13.91 ^^ gst 

n^iJl n (R X A)) > / •H°(J{' n (R X {t}))dt 

JteA 



> 



[ 2dt^2C\A). 
JteA 



On the other hand, the geometry of Jf implies that 

n\jt niRx A)) < V2C\A). 

As a consequence we get H^iJ^) > n^{Jf). The inequality n^{J^) > n^{Ji) 
follows in the same way. D 

We observe that dii{- ,951) is not in general an optimal solution of (Pi), as the 
following two examples show. 
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Example 3.10. Consider the subsets J7+ := {x e K^ : ||a;||ii < 3} and fl- := {x G 
M^ : ||a;- (2,2)||/i < 1} and set Q := mi(f7+ U fi-). Define u as 

u{x) =dii{x,dn+)xn_i_{x) - dii{x,dil^)xn.{x) ■ 

Then J"(w) = 16, while J"(d,i (• , dfl)) = 16 + 2^2. 

The fact that the function d;i(- ,dft) is not a solution of problem (Pi), is not 
a question of sign, as one could think from the previous example (notice that u is 
not positive on fl). Indeed in the next example we will exhibit a positive function 
V e £pwi^) such that J"(d;i (• , d^)) > T{v). 

Example 3.11. Define the four following squares Ci :— {x e K^ : ||a;||/i < 12}, 
C2 := {x e R2 : \\x - (10, 10)||,i < 8}, C3 := {x G M^ . ||^ „ (2, ll)||,i < 1}, 
Ci := {x G K^ . ii^. _ (11,2)11,1 < 1}. Let Q. := mt(Ci U C2 U C3 U C4). It is easy 
to see that J"(d,i (• , dVL)) = 120 + 16^2. Now, let us define 

u{x) = dii{x,d{Ci UC2))xciuC2(a;) + dii{x,dC3)xc3{x) + dll{x,^C4)xc4{x)■ 
Th.en u is non-negative, u is a solution of (jl.ip and -F(u) = 88 + 24-\/2, and thus 
dji(- , 9r2) is not optimal for (Pi) among the non-negative functions of £pw{^)- 

4. The general case 

We now turn to the general case where 51 is any bounded open connected subset 
of K''^ with Lipschitz boundary. In the previous section we saw that the functional 
J^ attains its minimum in £'pu,(51), if the set i7 satisfies assumption (H). In the 
general case we can't expect T to be finite in &{fl), as the following example 
shows. 

Example 4.1. Let 57 — (0,1)^: we claim that JF{v) = -fcx) for every v G &{ft). 
Indeed, let us consider, for t G (0, 1/2) fixed, the points (s,i) in ft with s G (0, 1). 
Then for almost every t G (0, 1/2) the function wt : s >—>■ v{s,t) has the following 
properties: w'^ = -§^{- -.t), \w[\ — 1 a.e., Wt{0) = wt{l) = and |w;t| < i, due to 
Lemma [3.41 This implies that w[ has at least [^] jumps, if [x] denote the integer 
part of a real number x. Therefore, applying Lemma 13.91 we have 



"1/2 j-1/2 

n\j^) > I n°{j^ n {{0,1) xt))dt > 



dt — -|-oo , 



andso J"(u) ^^T-C'{Jl 



In view of the last example, a possible way to generalize the analysis of the 
previous section is to measure the jump set of Du for u G &{^) through a "weight" 
function h G Co (57) which penalizes the jumps of Du in the interior of 57. 

We propose the following construction for the function h. Let (57„)„>i be a 
sequence of open subsets of 57 satisfying hypothesis {H) and approximating 57 in 
the following way: 

Vn>l ,\ , <rf(51„,951)<-. 

n + 1/2 n 
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We set rio = 0. Notice that the sequence (ri„)„>i is increasing, and that for any 
n > 1 the open set a;„ :— fin \ ^n-i satisfies hypothesis (H). Then for any n > 1, 
consider the problem 



It fohows from Theorem 13.11 that (5„ := inf(P") < +oo for any fixed n > 1. We 
now define h € Co(^, ]0, +oo[) by 

l/[4(Qi+a2)] if xe fii, 

dii (x, d^n) + dii (x, Orin+i) 



h{x) := < 



{n + l)2(a„ + a„+i) [n + 2)2(a„+i + a„+2)_ 



if a; e a;„+i, n > 1, 



where a„ := max{l, (5„ + H^ ^(9ri„)} for any n > 1. 
Define 

£(f7) := {v e e{n) : Dv E SBViodn)'^} ■ 

We win show the foUowing resuh: 

Theorem 4.2. Let fl be an open bounded connected subset of M.^ with Lipschitz 
boundary. Let h defined as above. Then the variational problem 

(Ph) inf I J'kiv) -^Y. I Hx)dn''-\x) : v G £{fl) I . 

has finite value and has a optimal solution. 

Remark 4.3. As for problem (Pi), (Ph) has at least two optimal solutions: if d is a 
solution, then —v is a solution too. 

Remark 4.4. If fl satisfies hypothesis (H), one can choose ft, = 1 as Theorem 13.11 
shows. 

Proof. We divide the proof into two steps: in the first one we show that inf (P^) is 
finite and in the second one we prove the existence of an optimal solution. 
Step 1: We show that there exists a function u € £"(17) such that J-h{u) is finite. 
For any n > 1 let ujn be the sets considered in the definition of h. Let u„ be an 
optimal solution of the associated problem (P"), extended by on R^ \a;„. (such a 
solution exists due to Theorem l3.ip . We then define u € £{n,) by u :— X^n^ Unix), 
and claim that J-h{u) < +oo. Indeed, using that h < — ^ on aJ„ for any n > 1, 
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one can estimate, for every i 



J? 



+00 

E 

n=l 



h{x)dn^-\x) 



J2'n(w„uaa„) 



+ OC 

E 



h{x)dH''-'{x) + 



h{x)dW'-'{x) 



j""nu 



Jrndn„ 



+00 


/ 5:'"'""'^-> 


n=l 


J""ncj,, 


+00 

> ^ < +00. 

n=l 



< 



< > -r 



J"n9n„ 



1 



dn'^-'{x) 



As a consequence inf(P;i) has finite value. 

Step 2: We show that there exists an optimal solution of problem {Ph ) applying the 
direct methods of the Calculus of Variations. Let (w")n C 8{Vl) be a minimizing 
sequence for {Ph)- As in the proof of Theorem 13. II we can assume that, up to a 
subsequence, u" — > v°° in Co{Vi) and weakly* in W^''^{Vi) for some v°° e Wo^'°°(r2). 
It remains to prove that v°° is an optimal solution of [Ph)- 

Let us prove that v°° € £{Vi). Since (w")n is a minimizing sequence we can 
assume that there exists C > such that J^h{v^) < C for every n € N. Fix an open 
subset uj oi D, with Lipschitz boundary and such that to G Cj G fl, and let a > be 
such that 0<a</i<linw. We can then write 



N 

^E 

1=1 






N-1 



(x), VneN. 



Let us fix i € {1, . . . , N}. Since ^- £ SBV{uj) and takes only the two values ±1 
for any n, we can then write 



c^u" 



dxi 



^ C\uj) + 2W^-\J^ nw) < C^<^) + 2C. 



BV{u 



We now apply Theorem 16.21 to the sequence {^^)n'- condition i) is the inequality 
given above; condition ii) is straightforward since V^^ = a.e. in w; for condition 
Hi) it is sufficient to choose / = 1, for example, because in this way 



Vne N, 



/ 



/" nw 



dv" 



dxi 



dU 



Af-l 



{x) 



n 



N-l 



(jfnLu) < a 



As a consequence S^ — )■ gi weakly* in BV{lu) for some function gi G SBV{uj) 



Since we already know that 



dxi 



— ^ 



dxi 



weakly* in L^{n), we infer that ^ 



dxi 



SBV{uj). Moreover ^ ^ ^ m L^{u}) so that |^| = 1 a.e. in w. This being 
true for any index i and any such open subset uj of J7, we infer that v°° belongs to 

£{n). 
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To prove that v°° is a minimizer of Jtj , let us show that 
(4.1) hminf / h{x)dn^-^{x)> [ h{x) dn^-\x) i ^ 1, . . . ,N 

n^oo J J 

for any io as above. To this aim we observe that for every n G N U {oo} one has 



h{x)d'H^-^{x) = ]- I h{x)d 



dxi 



(x). 



Therefore (|4.ip foUows from Theorem 16.41 This is sufficient to conclude the proof, 
since 

liminf f h{x)dn^-\x) > sup / h{x)dn^~\x) = / h{x) d'H'^-^x) 

where the supremum is taken over the open subsets uj d U d fl with Lipschitz 
boundary. We then infer from the last inequality that liminf J^/i(w") > J'h{v°"), 

that is, v°° is a minimizer of J^h- D 

5. Open problems 

We would like to present here some open problems related to our work: some 
questions about problems (Pi) and (P/j) and the question of the selection of a 
particular solution of the differential problem (|l.ip . 

1. Suppose that il satisfies property (H). We showed in section [3] that there exists 
a function u £ Cp^{i}) solution of problem (Pi). Can we say that the number of 
the connected components of Du is finite? This seems to be natural, as shown in 
Proposition 13.51 

Another natural question is: has problem (Pi) a finite number of solutions? Does 
it admit just two optimal solutions (a solution and its opposite)? This is true when 
f2 is a rectangle and TV = 2 as shown in Proposition 13.81 With the same arguments 
as in the proof of this proposition, it is easy to show that the functions ±m, where 
u is the solution defined in Example 13. 101 are the only optimal solutions of problem 

(Pi)- 

Another related question is to extend the result of Proposition 13.81 to higher 

dimension, that is to prove that if fi a convex open subset of R^ satisfying (H) 

then the distance functions ±d;i(., dH.) are the only solutions of problem (Pi). We 

notice that under the above hypothesis on f2, the function dii{.,dft) is a concave 

solution of ()1.1|) : this implies that on a non-void interval IR x ynil (where y G R^~^) 

the partial derivative gf-c'ii (•, 957) has only one jump. As a consequence the jump 

set of ■^rdii{.,dil) is also in some sense minimal since on each such interval it is 

reduced to one point, and so it has minimal cardinality. 

2. In sectionUwe studied the general case, that is when H. is any subset of R-^. A 
natural question is the study of the asymptotic behavior of a family {uh)h of the 
optimal solutions of (Ph) when /i — ^ 1 locally uniformly in 51. The limit problem is 
then obviously (Pi), which in general is ill-posed (since it has no feasible solution), 
but one expects uu to converge uniformly to a function u e £(il) having some local 
minimization properties in the open subsets w C cJ C fl. 
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3. In section |4] we saw that for a general open set il, J^(f) may be infinite for any 
V e &{fl). This means that for some index i one has H^^^{Jl') = oo. On the other 
hand we know that Ti,^ {J^) = 0, so a natural question is to identify the Hausdorff 
dimension s € [N — 1,N] for which H*{J^) is finite ii t < s, infinite if i > s. Then 
one may try to minimize this dimension s over ©(51), and then proceed as in section 
|3] with the optimal s in place of A^ — 1 . 

4. In this article we proposed to minimize the measure of the discontinuity sets 
of the gradients of the elements of ©(il) through the functionals J^ and J^h- Now, 
let {Gj}jg/ be the connected components of the Du~^{E) for u G ©(fl), where we 
recall that E = {x : Vi, \xi\ = — ^}- One could instead consider functionals of the 
form 

^5(H^-i(G.)) w^th g{t) = \ f' JJ| 

where /? e]0, 1[ . The advantage of this last type of functional with respect to our 
choice is that we can penalize the number of components Gi thanks to g. 

5. One of the aim of this work was to select a particular solution of (jl.l[) . The idea 
developed in the previous section is first to restrict ourselves to those solutions of 
(|l.ip that minimize problem (Ph)- Now how can we select a particular solution of 
{Ph) (up to its sign)? Does problem (Ph) admit more than two solutions? Notice 
that in any case the selection of these two solutions of our problem (jl.ip depends 
on the choice of h, so that it would be even more interesting to select a solution 
obtained by letting ft, go to 1 (see 4. above). 

In the case where problem (Ph) admits more than two solutions one can think 
of the following way to select a particular solution. Let Sq be the set of solutions 
of problem {Ph), and let (/„)„>i be an orthogonal basis of L^(ri) starting from 
fi — I. With the same arguments as in the proof of Theorem l4.2[ it is not difhcult 
to show that problem 



(Bi) sup<^ / ufi,ueS(] 



admits a solution. Let Si be the set of optimal solutions of problem {Bi ). li Si is 
composed by more then one function, we consider 



(-B2) sup<^ / uf2,ueSi 



and so on. If two functions u,v Cz Sq satisfy 



n Jn 

then one gets u = v. Therefore this method allows us to select a particular solution 
of (Ph). 

6. Appendix 

We recall here some compactness and semicontinuity results about BV, SBV 
and piecewise constant functions; for the proofs, see respectively Theorem 3.23 in 
[2], Theorem 1.4 in [1], Theorem 4.25 and Theorem 2.38 in ^. 
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Theorem 6.1 (compactness for BV functfons). Let Q, be a bounded open subset of 
M^ with Lipschitz boundary. Assume that {un)n is an uniformly bounded sequence 
in BV{fl). Then there exist a subsequence {un^)k o,nd a function u G BV{n) such 
that Urn. ^^ '^ weakly* in BV{Cl). 

Theorem 6.2 (compactness for SBV functfons). Let fl be a bounded open subset 

of M.^ with Lipschitz boundary. Let (un)n be a sequence of functions in SBV{i^) 

and assume that 

i) the functions u„ are uniformly bounded in BV{Q); 

ii) the gradients Vu„ are equi-integrable; 

Hi) there exists a function f : [0, oo) — >■ [0, oo] such that f{t)/t — > oo as t — > 0+ and 

f([ur,])dH^-^ < C < oo Vn e N. 

Then there exists a subsequence {un^)k and a function u G SBV{^) such that 
Urifc -^ u weakly* in BV{^), with the Lebesgue and jump parts of the derivatives 
converging separately, i.e. DaUn,. — >■ DaU and DjUn^^ — > DjU weakly* in Ai(n,K.^). 

Theorem 6.3 (compactness for piecewise constant functfons). Let fl be a bounded 
open subset of M.^ with Lipschitz boundary. Let («„)„ be a sequence of piecewise 
constant functions in il such that WunWi^ootm + H^^^iJun) is uniformly bounded. 
Then there exist a subsequence {un^)k o,nd a piecewise constant function u such 
that Un^ -^ u in measure. 

Theorem 6.4 (semfoontinuity in BV). Let il be a bounded open subset ofW^. Let 
iun)n be a sequence of functions in BV{fl) such that Un — > u weakly* in BV{fl). 
Then 

f{x)d\Dju\{x) < liminf / f{x)d\DjUn\{x) 
! "^°« Jn' 

for any non-negative continuous function / : O — > [0, +oo[ . 
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